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Abstract 

We study the well-known sociological phenomenon of gang aggregation and territory formation through an 
interacting agent system defined on a lattice. We introduce a two-gang Hamiltonian model where agents have 
red or blue afhliation but are otherwise indistinguishable. In this model, all interactions are indirect and occur 
only via graffiti markings, on-site as well as on nearest neighbor locations. We also allow for gang proliferation 
and graffiti suppression. Within the context of this model, we show that gang clustering and territory formation 
may arise under specific parameter choices and that a phase transition may occur between well-mixed, possibly 
dilute configurations and well separated, clustered ones. Using methods from statistical mechanics, we study 
the phase transition between these two qualitatively different scenarios. In the mean-fields rendition of this 
model, we identify parameter regimes where the transition is first or second order. In all cases, we have found 
that the transitions are a consequence solely of the gang to grafhti couplings, implying that direct gang to gang 
interactions are not strictly necessary for gang territory formation; in particular, graffiti may be the sole driving 
force behind gang clustering. We further discuss possible sociological - as well as ecological - ramifications of 
our results. 
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1. Introduction 

Lattice models have been extensively used in the physical sciences over the past decades to describe a wide 
variety of condensed matter equilibrium and non equilibrium phenomena see e.g., the reviews in [SI 1651 147] . 
Magnetization was the original application, but the list has grown to include structural transitions in DNA 
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[551 [TJ IM] , polymer coiling [SH [H] , cellular automata [3T1 [M] , and gene regulation [521 HI] to name a few. 
The resulting models are certainly simplified, but what they lack in detail is compensated by their amenability 
to analytical and computational treatment - and, occasionally, to exact solution. Moreover, at least for the 
behavior in the vicinity of a continuous transition, the simplifications inherent in these approximate models may 
be presumed to be inconsequential. In short, lattice models have proved extremely useful in the context of the 
physical, biological and even chemical sciences. In more recent years, lattice models have also been applied to 
study social phenomena [21 EHl HH] , such as racial segregation [53 EOj : voter preferences [32 EH US] , opinion 
formation in financial markets [59l [Ml 13], and language changes in society [37l |53l [56], offering insight into 
socioeconomic dynamics and equilibria. In this paper we consider the problem of gang aggregation via graffiti in 
what is - to the best of our knowledge - the first application of lattice model results to the emergence of gang 
territoriality. 

Scratching words or painting images on visible surfaces is certainly not a new phenomenon. Wall scribblings 
have survived from ancient times and have been used to reconstruct historical events and to understand societal 
attitudes and values. Today, graffiti (from the Italian graffiare, to scratch) is a pervasive characteristic of all 
metropolitan areas [5j. Several types of graffiti exist. Some are political in nature, expressing activist views 
against the current establishment; others are expressive or offensive manifestations on love, sex or race. At times, 
the graffiti is a mark of one's passage through a certain area, with prestige being attributed to the most prolific 
or creative tagger or to one who is able to reach inaccessible locations. The mark can be anything from a simple 
signature to a more elaborate decorative aerosol painting [SH |4] . All of these types of graffiti are usually scattered 
around the urban landscape and do not appear to follow any predetermined spatio-temporal pattern of evolution. 
They affect the quality of life simply as random defacement of property, although sometimes they are considered 
art [55] . 

On the other hand, gang graffiti represents a much more serious threat to the public, since it is usually a sign 
of the presence of criminal gangs engaged in illegal or underground activities such as drug trafficking or extortion 
[40[ I28j. Street gangs are extremely territorial, and aim to preserve economic interests and spheres of influence 
within the neighborhoods they control. A gang's "turf" is usually marked in a characteristic style, recognizable 
to members and antagonists [T51 135] with incursions by enemies often resulting in violent acts. The established 
boundaries between different gang factions are sometimes respected peacefully, but more often become contested 
locations where it is not uncommon for murders and assaults to occur [llj . It is here, on the boundaries between 
gang turfs, that the most intense graffiti activity is usually concentrated. 

Several criminological and geographical studies have been presented connecting gang graffiti and territoriality 
in American cities [5^1111115]. In particular, it is now considered well-established that the spatial extent of a gang's 
area of influence is strongly correlated to the spatial extent of that particular gang's graffiti style or language. 
Furthermore, it is known that the incidence of gang graffiti may change in time, reflecting specific occurrences 
or neighborhood changes. For example, rival gangs may alternate between periods of truce and hostility, the 
latter being triggered by arrests or shootings. Similarly, boundaries may shift locations when the racial or socio- 
economic makeup of a neighborhood changes, creating new tensions, or when gang members migrate to new 
communities [5]. In all these cases, periods of more intense gang hostility are usually accompanied by intense 
graffiti marking and erasing by rival factions in contested or newly settled boundary zones 15] . 

The purpose of this paper is to present a mathematical model that includes relevant sociological and geo- 
graphical information relating gang graffiti to gang activity. In particular, we study the segregation of individuals 
into well defined gang clusters as driven by gang graffiti, and the creation of boundaries between rival gangs. We 
use a spin system akin to a 2D lattice Ising model to formulate our problem through the language of statistical 
mechanics. In this context, the site variables Si have two constituents which represent 'gang' and 'graffiti' types, 
respectively, and phase separation is assumed to be the proxy for gang clustering. For the purpose of simplic- 
ity, we consider only two gangs, hereafter referred to as the red and blue gang, whose members we refer to as 
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agents. Lattice sites may be occupied by agents of either color or be void. Since gang members are assumed to 
tag their territory with graffiti of their same color, we also assign a graffiti index to each site representing the 
preponderance of red or blue markings. 

In particular, agents are attracted to sites with graffiti of their same color, and avoid locations marked by 
their opponents. We deliberately avoid including direct interactions between gang members, so that "ferromag- 
netic" type gang-gang attractions exist only insofar as they are mediated by the graffiti. On one hand this 
is mathematically interesting: in the broader context of physical systems, interactions are often mediated but 
rarely are indirect interactions the subject of mathematical analysis. On the other hand, by excluding direct 
gang interactions, we can specifically focus on the role of graffiti in gang dynamics and segregation. Furthermore, 
as will be later discussed, under certain conditions, gang-gang couplings may be unimportant, and one of the 
primary conclusions of this work is that they appear to be unnecessary to account for the observed phenomena 
of gang segregation. In any case, we informally state without proof that all the results of this work also hold if 
explicit agent-agent interactions are included. 

We thus write Si = {rji^gi), representing the agents and graffiti configuration at site i, respectively. The former 
component rji is discrete allowing, for simplicity, at most one agent on each site. The latter gi is continuous and, 
in principle, unbounded. We let s denote a spin configuration on the entire lattice, and in Section [2j propose 
a Hamiltonian, Jff{s), to embody all relevant sociological information. Once J^{s) has been determined, the 
probability for the occurrence of a spin configuration s on a finite connected lattice A C is determined by the 
corresponding Gibbs distribution F(s). Note that due to the choices made on the range of the ?7i, (?i values, F(s) 
is discrete in the rj variables and continuous in the g ones. It is given by 

F(s)-iexp(-,ir(s)), 
where Z is the partition function for the finite lattice A formally provided by the expression 

Z = ^exp(-jr(s)). 

Here, S denotes the set of all possible configurations on A and the summation symbol is understood to be a 
summation over the discrete components and an integration over the continuous ones. As usual, we begin with a 
finite lattice and its associated boundary conditions, and obtain infinite volume results by taking the appropriate 
limits. Using techniques from statistical mechanics, we prove that our system undergoes a phase transition as the 
coupling parameters are varied. In the unconstrained ensemble, certain parameter choices lead to predominance 
of either the red or blue gang, indicating that for configurations where the red to blue gang ratio is fixed at unity, 
a phase separation will occur. Conversely, in other regions of parameter space, there is no dominance of either 
gang type, indicating that the two are well-mixed and/or dilute. In this work we will investigate under which 
conditions to expect phase separation or gang dilution. 

Our paper is organized as follows: in Section |2j we give details of the model and in Section [3) we prove 
that a phase transition exists as a function of the relevant parameters. Since information on the location of all 
transition points is, by necessity, incomplete we consider an approximation in the form of a simplified mean field 
version of our Hamiltonian and derive the corresponding mean field equations in Section [4j Here, we show that 
the mean field Hamiltonian also exhibits a phase transition and we further prove that the latter is continuous in 
one specified region of parameter space and first order in another. Finally, in Section [5] we end with a discussion 
of potential sociological and ecological implications of our results. 
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2. The Hamiltonian 

Let us define a spin system on a finite lattice A C Z^. Here, the spin at each site i e A is denoted by 
Si = (jli,gi) and, we reiterate, rji denotes the agent spin and gi represents the graffiti field. We allow the agent 
spin to be in the set {0, ±1}; 77^ = — 1 if the agent at site i belongs to the blue gang, 77^ = if there is no agent, 
and rji = +1 if the agent is a red gang member. The graffiti field is in the set of real numbers: gi > indicates 
an excess of red graffiti, gi < an excess of blue graffiti, and, in either case, \gi\ indicates the magnitude of the 
excess. We now introduce the formal Hamiltonian ^(s) 

- Jf{s) ^JYI + "^^9^ + « I] - ^ II 5', (1) 

<i.j> i i i 

where s is a given configuration on the full A lattice, i and j index its sites and sum taken over 

every bond between nearest neighbor sites belonging to A. We discuss the role of spins on the lattice boundary 
A'^ in Proposition |3.4| and following sections. The expression in [T] will be referred to as the GI-Hamiltonian 
(graffiti interaction Hamiltonian) and its corresponding partition function will be denoted by an unadorned Z. 
Note that since rji is either or ±1, rjf — \rii\; however, we choose to display the above form to leave open the 
possibility of 77^ € Z. As discussed earlier, there are no explicit agent-agent interactions in this model; indeed, the 
structure of the Hamiltonian assumes that gang members interact with each other only via the graffiti tagging. 
As a result, occupation at site i by a gang member is "energetically" favored only if nearest-neighbor and on-site 
graffiti are predominantly of its same color. The two coupling constants, J for nearest-neighbor interactions and 
K for on-site occupation, reflect this trend. The arjf term represents the proclivity of a given site to be occupied 
by agents regardless of color, implying that gang members carry a strong tendency to occupy unclaimed turf if 
a ^ 1, while a <C — 1 represents a natural paucity of gangs altogether. Finally, we assume graffiti imbalance 
of either color to be energetically unfavorable via the —Xg^ term. This can be interpreted as natural decay of 
graffiti due to the elements, or to police or community intervention. For purposes of stability, A must be positive. 
Although the interactions J, K are tacitly assumed to be positive, generalizations to negative values may be 
possible, and a corresponding analysis may be undertaken given the proper sociological interpretations. 

3. Phase transition in the GI system 

3.1. Low temperature phase 

The basic strategy we follow to demonstrate an ordered, "low temperature" phase is a contour argument, here 
illustrated: Suppose that rji, the agent spin at site i, differs from the agent spin rjj at a different site j. The two 
agent spins can differ either by color, representing two different gang affiliations, or by occupation, where one 
site is occupied and the other is void. At the scale of nearest neighbors, each edge in the lattice can be defined 
as either a coherent bond, where the adjoining lattice sites are occupied and their agent spins are identical, or 
as an incoherent bond if this condition does not hold. Thus, explicitly, (rji, rjj) = (1, 1) or (—1, —1) are coherent, 
and all the other types are not. 

Let us now consider any path on the lattice that joins sites i and j. Since i and j have agent spins which are 
not identical, it must be the case that on any path between i and j, there is an incoherent bond. Furthermore, 
these incoherent bonds must form a closed contour on the dual lattice that separates i from j. In the following 
subsections, we derive a bound on the probability of any such incoherent bonds and their aggregation into 
contours. When these probabilities are small enough ~ which happens in certain regions of the parameter space 
- we can establish a low temperature phase. For example, the presence of a red agent at the origin will imply 
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that, with significant probabihty, the majority of the other sites will also be occupied by red agents, showing the 
existence of a red phase. Similarly, a blue phase can be shown to exist. 

To achieve all of these ends, we will employ the methods of reflection positivity described in [7] and |57j which 
contain a detailed account of useful techniques along with relevant classic references. In this paper, we will be 
working on the Lx L diagonal 2D torus - the SST - which we denote by T^. We will often refer to the Gibbsian 
probability measure on associated with the Hamiltonian in Eq.Q which we denote by Pl(-)- 

3.1.1. Reflection positivity 

By means of the reflection positivity of the Gibbs distribution we can easily bound the expectation of an 
observable which depends only on the spin at any two neighboring lattice points. This result will be used to 
build the contour argument that will lead us to prove the existence of a low temperature phase. We thus briefly 
introduce the concept of reflection positivity, referring the interested reader to [7] for a more detailed discussion 
of these topics. 

Consider a plane of reflection p which intersects the torus in a path running through next nearest (diagonal) 
pairs of sites. Let dp be the reflection operator through p. On the SST, this plane p divides the lattice into two 
halves, identified as and T^, such that n = p. Let denote the set of functions which depend only 
on the spin variables in and similarly for '^p~ . The reflection map, dp, which, in a natural fashion identifies 
sites in with those in T£ via a reflection through p, can also be used to define maps between '^^^ and '^p~: 
Specifically, if / € "^p^, we define dpf e '^p' to be the function / evaluated on the configuration reflected from 

A measure /i is reflection positive with respect to dp if for every f,g(z "^p , or , the following two 
properties hold 

1. E,,(/i9p/) > 0, 

2. E^(/V)=E^Mp/)- 

It is known (e.g., see [7]) that is reflection positive with respect to dp for every p of the above described type. 
We next use reflection positivity to find an upper bound on the expectation of observables defined on bonds. In 
doing so, we use the following lemmas: 

Lemma 3.1 Let {i,j) denote a bond ofTj^ and let and denote site events at the respective endpoints of 
the bond. Let zi^^'^^ denote the partition function ( on ) which has been constrained so that at each site with 
the parity of i, the translation of the event ai occurs and similarly for 7. Then, for i = 2*^ for some integer k, 



< 



r(a.7)' 



where V = L^ is the volume of the torus. 

Proof. The result from this Lemma dates back to the original papers on the subject. In particular, the use of 
bond events on the SST was highlighted in A modern and complete derivation is contained in j^. Section 
5.3. □ 

For a slightly more general scenario, let us consider the bond {i,j) and various events al,jj, . . . , a", 7" and let 
us denote by 61 — aj (17^ ... 5„ = a" n 7" the corresponding bond events as described. Letting b — U"^]^6j we 
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find 



Finally, we have 



r("3.7j) 



2V „ 



Lemma 3.2 Lei ri, . . . denote translations of the bond and br- the translation of the bond event(s) b 

described above. Then 



.(nr=A,)< 



Proof. Again, we refer the reader to [7], Section 5.3. 



E 



□ 



3.1.2. A bound on the incoherent bond probabilities 

In order to prove a phase transition by a contour argument, we must place an upper bound on the probability 
for the occurrence of any type of incoherent bond where agent spins of neighboring sites are different. There are 
four types of incoherent bonds, namely {r]i,r]j) = (— 1, 1), (— 1, 0), and (1,0), and (0,0), regardless of order. Let 
us introduce the following notation: consider undirected bonds between two particular neighboring lattice sites, 
{i,j) and let (•, •) denote the event of any of the nine coherent or incoherent bonds so that 



(•, •) e {(+, +), (-, -), (+, -), (-, +), (+, 0), (0, +), (-, 0), (0, -), (0, 0)}. 



Similarly, let Zj'^ ^ denote the partition function restricted to configurations where all agent spins are frozen in 
accord with the above described (chessboard) pattern and the rest of the statistical mechanics is provided by the 
graffiti field against this background [? 32J. The following is readily obtained: 

Proposition 3.3 The above described (agent-constrained) partition functions are given by 



A-,-) 



A+.-) 



= Z. 



r(0,0) 



(+■+) 



e \/7r 1 



[4J+Kf 



62 VTI" J_K^ 



-esA 



Proof. Since the agent variables are frozen, the gi Gaussian variables are independent and the above amount to 
straightforward Gaussian integrations. □ 



Using Lemma 



3.1 



and the fact that the full partition function satisfies Zj^ > Z, 



(+,+) 



we can write 
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Pl(0,0) < e-5"e-s^[4'^+^l', 
Pl(+,-)-Pl(-,+) < c-^, (2) 
Pl(+,0)----=Pl(0,-) < e-i"e ^ . 

We denote by e = e( J, K, A, a) the sum of the estimates for the probabilities provided by the right hand sides 
of the preceding display. For fixed a and > 0, note that as J\~^^^ — ?► oo (or, better yet, J\~^^^ and KX^^^^ 
both tending to infinity) the quantity e tends to zero. This implies the suppression of all incoherent bonds so that 
the lattice must be almost fully tiled with coherent ones. In particular, the lattice is nearly filled with agents, 
which, at least locally, are mostly of the same type. As will be demonstrated below, this implies the existence of 
distinctive red and blue phases, i.e., in the language of statistical mechanics, of a "low temperature" regime. We 
formalize this result in the next subsection. 

3.1.3. The contour argument 

We have now established all the tools we need to complete the contour argument. Accordingly, we now show 
that two well-separated lattice sites must, with probability tending to one, have identical agent spins in the limit 
e <C 1. This in turn will imply the existence of a low temperature phase. 

Theorem 3.1. Consider the Gl-system on 1? and let e{J,K,\,a) denote the quantity described in the last 
paragraph of the previous subsection. Then, if the parameters are such that e is sufficiently small, there are 
at least two distinct limiting Gibbs states characterized, respectively, by the abundance of red agents and the 
abundance of blue agents. Moreover, this property holds in any limiting shift invariant Gibbs state. 

Proof. Let us start on T/, with L = 2^ . For i,j e where i and j are well separated, let us consider the event 
Vb '■= {'Hi 7^ rjj} U {rji = 0}. We will show, under the stated conditions, that uniformly in L this probability 
vanishes as e -> 0. As discussed previously, in order for this event to occur, the sites i and j must be separated 
by a closed contour consisting of bonds dual to incoherent bonds. For £ — 4,6, . . . let = ^iii — j, L) denote 
the number of such contours of length ^ on T^. Then we claim that uniformly in L and i — j, 

me < 2fXi 

where A2 (with A2 ~ 2.638... < 3) is the connectivity constant for [35]. A word of explanation may be in 
order. The A2 generously accounts for walks of length £ in the vicinity of site i and the factor of two for walks in 
the vicinity of site j. Finally, the factor of i'^ accounts for the origin of the walk. Note this is an over-counting, 
e.g., contours which wind the torus but do not necessarily "enclose" i or j are counted twice. Using Lemma |3.2| 
we may now write 

The above obviously tends to zero as e — ^ demonstrating that in finite volume, the lattice is either populated 
with mostly red agents or mostly blue agents depending - with high probability - on what is seen at the origin. 
The implication of this result is that, for e sufficiently small, there are at least two infinite volume Gibbs states 
- which can be realized as the limits of the appropriately conditioned T^'s. These states have one of the two 
mutually exclusive characteristics: a preponderance of red agents or a preponderance of blue agents. The fact 
that the above must also hold in any shift-invariant Gibbs state is the subject of Theorem 2.5 and its Corollary 
in [5j with a slight extension provided by Corollary 5.8 in [T3]. □ 
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3.2. High temperature phase 

As is sometimes (e.g., historically) the case in statistical mechanics, it can be an intricate job to establish a 
high temperature phase - a region of parameters where the limiting Gibbs measure is unique and correlations 
decay rapidly. Typically, one calls upon the Dobrushin uniqueness criterion |23| . However for us, this route is 
interdicted by the unbounded nature of the gi graffiti field. The strategy here will be percolation based: First 
we establish the so-called FKG property for all the associated Gibbs measures. Then what follows will be a 
relatively standard argument through which we show that the necessary and sufficient condition for uniqueness 
is that the average of rji - akin to a magnetization - vanishes in the state designed to optimize this quantity. 
Then, finally, we will develop a random cluster-type expansion demonstrating that under the expected high- 
temperature conditions for the couplings, e.g., A 3> f , the stated condition on this magnetization is satisfied. In 
addition, high-temperature behavior should also be achieved under the condition that agents are sparse. This 
requires an alternative percolation criterion used in conjunction with the above mentioned expansion. In both 
scenarios, the rapid decay of correlations arises as an automatic byproduct. 

3.2.1. FKG properties 

In this paragraph we will demonstrate that the FKG Lattice Gondition (see e.g., [H] Page 78) is satisfied by 
any finite volume Gibbs measure associated with the GI-Hamiltonian. Let us start by noting that we can define 
a natural partial ordering on the pair of states Sj and s[ via the notation 

Si h si if rii > 77- and .g,, > 5-. 

Further we introduce the notation s ^ s' to signify that the above holds for all the Si, s'^ at each i G A. For 
individual spins Si and s^, we also denote Si V := (maxjyy^, 77^}, maxj^^, g^}) and similarly for the "minimum" 
Si A S-. Finally, for spin configurations s and s', the configurations s V s' and s A s' are defined as the sitewise 
maximum and minimum, respectively. The FKG lattice condition - conveniently stated for finite volume measures 
- is that for all s, s', the following inequality holds: 

Fa(s V s')Fa(s a s') > Fa(s)Fa(s'). (3) 

The well known consequence of the above is that any pair of random variables that are both increasing with 
respect to the partial order described above are positively correlated. 

Proposition 3.4 The finite volume Gibbs measures associated with the GI-Hamiltonian satisfy the FKG lattice 
condition. 

Proof. We consider an arbitrary graph and, as will be made evident, the proof automatically accounts for any 
fixed boundary conditions. Now, as is well known, it is sufficient to establish that the lattice condition Eq.(|3| 
holds when differences between configurations are exhibited only on a pair of spin- variables. The fixed boundary 
spins thus may be regarded as part of the background which is common to all four possible agent-graffiti spin 
configurations in question. Let us thus assume that the differences between two configurations occur at sites 
a and b in the graph where certain specified variables have been "raised" above a base configuration level s. 
We denote the single raise configurations by Sa and s;, and the double raise by Sab- Thus, it is sufficient to 
show F(sa6)F(s) > F(sa)F(sfc). All told, there are three possibilities to consider: graffiti-graffiti, gang-graffiti 
and gang-gang raises on the a and b sites. For the mixed gang-graffiti case we must also consider the a = b 
possibility where the gang and graffiti spins have been "raised" at the same site. We need not consider the 
normalization constant in any of these cases, since it appears in identical roles on both sides of the purported 
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inequality; consideration of the Boltzmann factors is sufficient. Let us introduce, in the setting of our general 
graph, the interaction 

- jr(s) = ^ Ji,jVi9j - ^[atvf + K9i], 
{hi) » 

where the first sum now extends over all edges considered to be part of the graph and our only stipulation is 
that Jij > 0. Also, we may formally include « = j in this sum. Let us denote the "raised" grafhti variables 
via the positive increments Sga and Sgb so that, in the graffiti-graffiti case, at sites a and b Qa ^ Qa + Sga and 
.95 9b + Sgij. It is straightforward to see that ,j^{sab) + --^(s) = ,3^{sa) + ^(sb) and the desired inequality 
holds as an identity. Similarly for the gang-gang case. We can now consider the mixed case where, without loss 
of generality, ga ^ ga + Sga and r/b ^ r]b + 5r]b, and for us, 5r}b = 1. Here 

-(=^(Sa) - ^(s)) = ^ Ji,a11iS9a + Ja,bmSga - >^5ga^ , 
ii^b 

while 

-(J^(s6) - ^(s)) = ^ JbjgjSrjb + Ja,bgJvb + aSrjb'^. 

However 

- {-^{Sab) - =^(s)) = J^,aViS9a + Ja,bVb5ga " AJfila^] + JbjgjSVb + Ja,b9aSVb + + Ja,bSga5Vb 

i^b jjtb 

> 2Jf (s) - JfiSa) - Jif{Sb). 

Combining the above results we find that indeed 

Jif{Sab) + ^(S) < Jf{Sa) + jr(S6). 

The same inequality can be easily shown in the mixed gang-grafhti case for a = 6, by assuming 5a ^ ffa + ^ffo 
and r]a rja + Srja and by following the same steps as above. This completes the proof. □ 

As an immediate consequence, we can identify boundary conditions on A which most favor the dominance of the 
red gang. Indeed, it is now seen - as was anyway clear heuristically - that we must make the boundary spins "as 
red as possible" in order for a predominance of A sites to be occupied by red agents. This amounts, somewhat 
informally, to setting gi = +oo and rji = 1 (which is anyway automatic if isT ^ 0) all along the boundary. This 
"specification" which seems a bit arduous to work with is not nearly as drastic as it sounds. Let us start with 
some notation: For A a finite subset of Z^, let us define dA as those sites in A'' with a neighbor in A and dA 
as those sites in A with a neighbor in A''. Clearly the only immediate consequence of the "drastic" boundary 
condition is to force = 1 for i e dA and to bias, by at most {3J + K)gi the a priori Gaussian distribution of the 
g's. We shall do that - and a bit more on dA arguing that this, at most, is the result of the "drastic" boundary 
condition on d{A U dA). Precisely, we define the red boundary condition on A as Tji = 1 and gi independently 
distributed as normal random variables with variance 1/2 A and mean (4 J + K)/2X for each i e dA. By the 
established monotonicity properties these are exactly the boundary conditions imposed on the slightly larger 
lattice that will optimize the average of rji and g^ for any i € A. 
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3.2.2. A uniqueness criterion 

It is not hard to show, by monotonicity, that a limiting red measure exists along any thermodynamic sequence of 
volumes and that the limit is independent of the sequence and therefore translation invariant. We shall denote 
this measure by /iR(-) and by Er(-) the corresponding expectations. Similarly for the blue measure we introduce 
^b{ ) and Eb(')- We can thus state 

Proposition 3.5 The necessary and sufficient condition for uniqueness among the limiting Gibbs states for the 
Gl-system is that Kn^rjo) = 0, where tjq is the spin at the lattice origin. 

Proof. For two measures and ^2 e.g., on {—1,0, 1}^ , we use the notation /ii > /i2 to indicate that for any 
random variable X which is increasing in all coordinates, the expected values Ei(X), calculated via the ^1 
measure are always greater than the those obtained via 112'- 

Ei(X) >E2(X). 

This is known as stochastic dominance. Consider /iR,(-) which, by slight abuse of notation, we temporarily take 
to be the restriction of to agent events. Suppose that Er(77o) — 0. Then, by translation invariance, we 
have V,^^{r]i) = for all i. Similar considerations apply to the corresponding /iB(')- It is immediately clear - by 
symmetry or stochastic dominance - that Pr(77o = 0) = Pb(?7o = 0) and thus the single site distributions are 
identical. By the corollary to the Strassen theorem [601116] since /zr > and these measures have identical single 
site distributions they must be identical probability measures. Similar considerations apply to the full measures 
since the distribution of the gi is determined by their conditional distributions given the local configuration of 
the Tyi's. Uniqueness is established since, if ij,q{-) denotes any other infinite volume measure associated to the 
GI-Hamiltonian, we have ^r > fiQ > /xb which implies equality in light of /ir — (J-b- D 

3.2.3. Proof of a high-temperature phase 

We shall develop a graphical representation for the Gl-system akin to the FK representation for the Potts 
model [29] that, for all intents and purposes, is the same as the one used in [10], where only the case of bounded 
fields is explicitly analyzed. Let us then consider the GI-Hamiltonian in finite volume with all notation pertaining 
to boundary conditions temporarily suppressed. For fixed s, we may decompose the graffiti fields and agents 
according to affiliation: 

9i = Qi-^f, i^i = ±1, Qi = \9i\, 
rji = riai, ai = ±1, = \rii\, 

where the ct's and i?'s have the definitive character of Ising variables. We can now write 

e^..9.'7. =e-^--«''-^(i?,,,,5^^,,^ +1), 

where Rij = R{Jij,qi,rj) := e^-'i^in^j — 1. In our case, we have Jij = J if i and j are neighboring pairs and 
Ji^i = K] which we will not yet distinguish notationally and consider a general Jj j label. Thus 

g-^(s) ^ -Q e-^'-9''-^ (i?,,,5^,,,^. + 1). 

Opening the product, we select one term for each "edge": If the Rij term is selected, we declare the edge to 
be occupied, otherwise it is vacant. It is noted here that the edges should be interpreted as directed: all edges 
appear twice and we must regard {i,j) as distinctive from {j,i); moreover, for K ^ 0, the above is understood to 
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include i — j. The configurations of occupied edges will, generically, be denoted by uj. Summing over the Ising 
variables, we acquire the weights 

where, as before the summation notation also indicates integration over the continuous variables. In the above, 
C{u)) denotes the number of connected components of uj; here connectivity deduced according to the directed 
nature of the edges or via a double covering of the lattice. Normalizing these weights by the partition function 
we obtain a probability measure on the bond configurations uj. As will be made explicit below, this probability 
measure on bond configurations is well defined in finite volume. Let us denote the probability measure on the 
bond configurations uj by P^(-), where the now denotes boundary conditions accounted for in a routine fashion. 
Then, for each uj consisting of appropriate edges, P^(i^) G (0, 1). We shall not discuss the problem of infinite 
volume limits which would take us too far astray but be content with statements that are uniform in volume. 
With regards to the latter, and of crucial importance for our purposes is the connection back to the spin-measure 
inherent in this representation. For the Potts models, this was first elucidated in [2 with the complete picture 
emerging in 25,. In particular, for any site, the contribution to the magnetization vanishes if the site belongs 
to a cluster that is isolated from the boundary. The principal objective for this representation is the following 
claim: 

Proposition 3.6 Let A C be a finite connected set and consider the above described representation in A with 
boundary condition on dA. Let (a, b) denote an edge with a ^ b and both a and b not belonging to dA. Let eat 
denote the event that this edge is occupied and let ui denote a configuration on the compliment of (a, b) . Then, 
for fixed a and K, there is an e(J, A) with £ — >■ as J^/A — > such that uniformly in A, uj and Q - as well as 
K and a, 

P?(e..) < 

Proof. Let W^(-) denote the configurational weights with associated boundary conditions as described above. 
Then it is seen that 

pQ(cjVeab) ^ VKQ(cjVeafc) 
l-P®(wVe,b) Wfiuj) 
Our goal is to estimate the right hand side of the above which thereby generates the quantity e featured in the 
statement of this proposition. Noting the positivity and product structure of the numerator and denominator, 
we may regard the object on the right as the expectation with respect to a weighted measure of the quantity Rah 
and we shall denote this by ]E^(i?Q^h). The latter will be estimated via conditional expectation: Let Qq^ denote 
a specification of the g-fields and agent occupation variables except for and let 

e := sup E„(i?afc I QgJ- 

Obviously, e > sup^ Ei^(i?Qb). As for the complimentary fields, there is not a great deal of dependence: In 
particular, all that is needed is that Tc = 1 for all c such that (a, c) € uj. Concerning the optimizing w, non-local 
considerations dictate simply, that uj be such that (a, b) does not reduce the number of components. Locally, 
as can be explicitly checked, or derived from monotonicity principals, the optimal scenario is when all bonds 
emanating from a are present in the configuration. Thus we have 

_ / e-^^-^+^'>''R*{J)R{K)e-^i^ dq _ J e"^?" sinh^( Jg) smh.{Kq)dq 
^ ~ J e-(3-'+i<')''R\J)R{K)e-^i' dq ~ J e-^i' smh"^ {Jq) smh{Kq)dq' 
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where in the first line R{J) := R{J,q, 1). We claim that the final ratio is bounded by J/A^/^ multiplied by a 
constant that may be proportional to the ratio K/X^^'^. Indeed let us substitute w = X^^^q and k := K/X^^^. 
The above quantity can thus be rewritten as 



e " sinhKw(sinha;J/A)^(iw 
J e~"^ sinh Kw(sinha; J/A)'^(iw 

Our claim is obvious if k — > but we may wish to consider cases where n stays bounded away from zero. In 
general, the integrands are not dominated by large lj and we may expand the factors sinhwJ/A with the result 

. 2J /e 



^1/2 J g-w2 gjnh^^ . ^3flijj 

We finally claim is that the right side is bounded by a linear function of k: 

1 /c- " 



1 + K J c sinh KLu ■ Lo^duj 



<B, 



for some B < oo. This is indeed true as k — >■ 0. We only need to show that the inequality holds in the case 
K — > oo. But here the factor e~" sinhKo; is, essentially, a Gaussian in the variable uj — k and the desired result 
follows. □ 



Theorem 3.2. Consider the Gl-system and let e denote the quantity described in Proposition 3.6 Then for 
e < So, given by 



1 2 1 



there is a unique limiting Gibbs state featuring rapid decay of correlations. 



Proof. Using the result of Proposition 3.6 we shall compare the described graphical representation with inde- 
pendent bond percolation on Z^. We start with a well known - and readily derivable - result: Let Yi, . . . Yat 
denote an array of Bernoulli random variables with collective behavior described by the measure /iy and let Yy 
denote a configuration on the compliment of Yj . Let us now introduce 

= maxPY(y,- 1 I Yy.) 

to denote the maximal conditional probability of observing {Yj — 1}. Finally, let Xi,. . .Xjy denote a collection 
of independent Bernoulli random variables with parameters pi , . . . ,Pn- Then, denoting the independent measure 
by MX, we have 

MX > My- 

Thus we may bound the probabilities of increasing events in the graphical representation by the corresponding 
probabilities from independent percolation on with bond occupation probabilities determined by the e from 
Proposition |3.6[ However, we must note that the relevant percolation problem has multiple types of edges. The 
Eq in the statement of this proposition bounds the probability of the event eat U ef,a by ^. If the clusters of 
the featured representation fail to percolate, then, as A Z^, the origin is disconnected from the boundary 



with a probability tending to one. As discussed just prior to Proposition 3.6 this implies Er (7/o) = and by 
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Proposition 3.5 uniqueness is established. Under the condition e < eo, exponential decay of correlations can also 
be established. We will be content with the decay of the two point function. The problem of general correlations 
under these conditions has been treated elsewhere [HI [17]. In particular, for i,j G Z^, £(77^77^) in the unique 
infinite volume measure is bounded, in finite volume approximations by the probability that i and j reside in the 
same cluster. For e < Eq^ this decays exponentially in \i — j\ uniformly in A for |A| sufficiently large. □ 

We now turn our attention to an alternative criterion for high temperature behavior which may also be of 
relevance in a sociological context; Sparsity of agents. Mathematically, this pertains to the situation where a is 
large and negative {—a 3> 1) which a priori suppresses the fraction of agent occupied sites. Our arguments will 
initially be based on more primitive notions of percolation and, following the methods of [20] (see also [T8l \19\ ) 
could, perhaps, be completed along these lines. However, it turns out to be far simpler to appeal to the graphical 
representation just employed for the final stage of the argument. We start with the relevant notion of percolation 
and connection. In the context of site percolation on Z^, we may define various notions of connectivity |34j . 
Here we define o-connectivity to indicate connection between sites that are no more than two lattice sites away. 
This is not to be confused with *-connectivity which does not consider a pair of sites to be connected if they are 
separated by two units in the vertical or horizontal direction. We denote by the threshold for o-percolation on 
Z^. Standard arguments dating to the beginning of the subject show that G (0, 1); in particular, p'^ is less than 
the threshold for ordinary, or even ^-connected, percolation and mean-field type bounds readily demonstrate 
that > j^. 

The next proposition concerns the relative abundance of, e.g., red sites under the condition —a 3> 1 with the 
other parameters fixed. 

Proposition 3.7 Consider the Gl-system with parameters \, K and J fixed. Then there is a — Sa{ J, K, X) 
with 6a ^ as a ^ —00 such that uniformly in volume and boundary conditions, for any site i that is away 
from the boundary 

P®(r;, = 1) < 6a. 



Proof. Here we employ the preliminary (red ^ blue) FKG properties that were established earlier, in 3.2.1 We 
start with a 7 > (and somewhat "large") and, for j G A not too near the boundary, we consider P'^igj > 7). 
By the FKG property, this probability is less than the corresponding conditional one given that rjj — 1 and that 
r]k — 1 for all k that are neighbors of j. This conditional probability is given by a definitive expression: 

The above can be expressed directly via the error function but in any case, as is not hard to show, 

"7—2 ' 

as long as 7 > (4J + K)/2X, which also quantifies how large 7 must be. Provided i is a few spaces away from 
the boundary, we note that 1 — 56-y is a valid estimate of the probability that both gi and the 5-values at the 
neighbors of i do not exceed 7. Let us denote this (good) non-high field event by Gi. Then we may write 

P®(r7, = 1) = P®(G,)P?(r?, = 1 | G.) + P® (G^P^t?, = 1 | G^) 
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where, in various stages we have employed worst case scenarios. Clearly, for fixed (J, K, A) we may choose 7 large 
so that 6^ is small, and a negative and large in magnitude, so that 5a is small. □ 

Theorem 3.3. Consider the Gl-System and suppose that —a is large enough so that 5a < pt described just 
prior to the statement of Proposition |3.7| Then there is a unique limiting Gibbs state featuring rapid decay of 
correlations. 



Proof. By the dominance principle stated at the beginning of the proof of Proposition 3.2 if < p*, the red 
agents fail to exhibit o-percolation regardless of boundary conditions. Now consider, in the context of the bond- 
representation, the event that the origin is connected to dA in the red boundary conditions, which represents 
the sole non-vanishing contribution to E^{rjo = +1). The bonds of any path connecting the origin to dA within 
this cluster may be envisioned as alternating connections between agents and fields; the connection to the red 
boundary ensures that both types of entities take on the red color. In particular, all the agents in the cluster are 
red so that these agents must (at least) form a o-connected cluster. Hence, in finite volume, we may bound 

Kivo = +1) < p;^(o dA) 

where {0 OA} is the event of a red o-connection between the origin and the boundary. 

o,R 

When the red agent occupation probabilities are dominated by independent sites with parameter 5a < P*' 
such probabilities decay exponentially. Evidently, in the limiting state, the "magnetization" vanishes which by 
Proposition |3.5| implies a unique state. Similarly, exponential decay of correlations is implied by exponential 
decay of o-connectivities. □ 



4. The Mean Field Rendition 

In the previous section, we showed that a phase transition between well-mixed and clustering configurations exists 
for the general Hamiltonian in Eq. ([T]). However, finding the exact or even approximate values of the J, K, a, A 
parameters for which the well-mixed to clustering transition occurs is in general a difficult task. Moreover, the 
nature of the transition is not elucidated by the techniques of the preceding section. On the basis of informal 
simulations described in the Appendix and certain other considerations it appears that the transition may be 
discontinuous or second order depending on where the phase boundary is crossed. This cannot be proved in the 
context of the present model. We thus introduce a mean-field Hamiltonian, where instead of nearest-neighbor 
interactions we consider an all-to-all (interaction) coupling that is rescaled by the number of sites. Models of 
this sort are often referred to as complete- graph systems. The mean-field Hamiltonian allows us to define, in the 
thermodynamic limit, a simple mean field free energy per particle. This free energy can be subjected to exact 
mathematical analysis which provides a quantification of the phase transition. In particular, we have found that 
the phase boundary between the diffuse states and the gang-symmetry broken phase can indeed be of either 
type. 

Let us thus consider a lattice of N sites - where the detailed geometry is no longer of relevance. At each site 
i, there is the same Si = (rii,gi) featured in the previous section. However now, the Hamiltonian reads 

- ^^^''^(s) = ^ E J^^9, + J2(^vf - ^gf). (5) 

It is observed that the couplings J and K need no longer be distinguished. Indeed, for large N , the girji interaction, 
and any other particular interaction is not of pertinence. We now introduce the relevant collective quantities, n 
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G, and b, obtained via rji and g^, which will allow for a more convenient analysis. In particular, if A''+ and N 
designate the number of red and blue lattice agents, respectively, we define 

N+ + N- , N+-N- 

:= — and n := 



N N 

as the fraction of the lattice covered by agents of any type and the excess - positive or negative - of this fraction 
that is of the red type. Moreover, wc introduce G = jj gi to be the average graffiti imbalance. In this 
context, n and G are akin to magnetizations in a standard one-component spin model, with n corresponding to 
magnetization in the agent variables and G in the graffiti field. For occasional use, we also define = N^/N < 1. 
We remark that in these definitions there is an implicit A'' dependence which is notationally suppressed. 

4-1- The partition function 

In the forthcoming, we will evaluate, asymptotically, the mean field partition function defined in accord 

with the previous section as the partition sum Zj^^ = e~'^^ Here, for reasons which will soon become 
clear, we will treat the graffiti field variables slightly differently. We define 

as the normalized Gaussian measure for the individual field variables. Letting g denote the array of these random 
variables we may write 



crMF Iff 



where Eg(-) denotes expectation with respect to the free (independent) ensemble of Gaussian random variables 
and denotes the rest of the partition sum i.e., over the agent configurations. It is acknowledged that this 
differs from the prior definitions by a multiplicative factor of [A/tt]^/^ which, of course, is inconsequential. 

It is at this point, with the current formulation, that the advantage of the all-to-all coupling is manifest: 
For any s (and any A^) the quantity in the exponent depends only on n, h and G: Z^^ = Eg(X]r; e^^'^"''-''^"'^^). 
Concerning the agent configurations, to perform the summation, we must multiply the integrand by the number 
of ways of arranging 7V+ red sites and A''" blue sites among A'' possible positions. We denote this object by 
WN{b,n) which is given, explicitly, by the trinomial factor 



As for the graffiti field configurations, it is noted that since G is proportional to a sum of Gaussian random 

variables, it is itself a Gaussian. Indeed the mean of NG is zero and the variance is A/'[2A]~"^. Thus the 
expectation over g can be replaced with the expectation over NG leading to 

= ^EjvG[W^iv(&,n)e^[''''G+«''l] oc ^ W^;v(n, 6)e^[''"^+«^-^^'ldG 

n,b n,b 

with the constant of proportionaUty independent of N. Now, on the basis of the Stirhng approximation, 



WN{b,n) 



b±n ^ -I -N 

b + n \ 2 f h — 



2 



(1-6) 



i-f) 
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Thus, modulo lower order terms, we have Zj^^ « ^„ ^ ^ e ^^(f>^'^,G) ^Y^^^e $, the free energy function, is given 

by 



b + n 



b-' 



(l-b) 



1-6 



(6) 



In accordance with standard asymptotic analysis 



lim ^ 

N—>-oo 



1 



logZj^^^ = min $(6, n, G) := Fmf 

I\ b,n,G 



where Fmf = Fmf{J, ct, X) is the (actual) limiting free energy per site. While various aspects of the above 
scenario for all~to-all coupling models have been long known and certain cases explicitly proven |26) . there is a 
general theorem to this effect that is sufficient for our purposes, presented in Section 5 of 9]. Thus the efforts of 
a mean-field analysis may be summarized as follows: we are to minimize $(6, n, G) and the values of b, n and G 
at the minima - as a function of the couplings - will determine the various phases of the system. Even in this 
simplified context, as will be seen, the phase transitions can be dramatic. 



4-2. The mean-field equations 

The free energy function is obviously well behaved except at the extreme values of the variables. In particular, 
we would like to assume that < & < 1 and —b<n< +b where the strict inequalities imply that the function is 
smooth. Now a direct calculation of the asymptotics makes it clear that no minimum could possibly occur near 
the 6 = 0, 1 and 6 = ±n boundaries. Thus, we can confine attention to the interior of the above b and n intervals 
and proceed by differentiation of ^{b, n, G) as defined in Eq.(|6]). Thus we arrive at the mean-field equations: 

9$ 



dG 
"db 
dn 



Jn - 2XG = 0, 



1 



a + log(l -b)- - log 



62- 



= 0, 



= JG 



1 



log 



= 0, 



(7) 
(8) 
(9) 



Free energy minimization only occurs for values of (n, b, G) that satisfy the above system. However, other 
stationary points for $ can - and, e.g., in the case of discontinuous transitions generically will - occur so we 



we rewrite Eqs. ([8| and ^ as 



must proceed with some caution. It is noted that Eq.(|7|) allows us to eliminate G altogether. Defining M = 



4e 



2a 



62 



fin 



(1-6)2 



- log 
2 ^ 



6-' 



(10) 
(11) 



The analysis of this system, along with the minimization it is supposed to imply will constitute the bulk of the 
remainder of this work. Foremost, it is noted that the presentation in Eqs. (10 1 and (11) arc, for all intents and 



purposes, the same as would have been obtained from the mean-field version of the so-called BEG model [T^]. 
As such, some aspects of the current problem have been treated in [57]. However, the specifics in [23 are not 
readily translated into that of the current work and, moreover, our conclusions are achieved by straightforward 
methods of analysis. 
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Our investigation will proceed as follows: It is evident from physical considerations, and the subject of an 
elementary mathematical theorem proved at the end of this subsection, that as the parameters sweep through 
their allowed values, a phase transition occurs from the circumstances where $ is minimized by n = to those 
where n ^ is required. First, we will follow the consequences of the assumption that this happens continuously: 
i.e., that the minimizing n goes to zero continuously through small values. In the leading order, this provides 
a purported phase boundary which we denote by the LSP-curve. Considerations of higher order terms in the 
vicinity of the LSP-curve yield that for certain portions of the curve, the stipulation is self-consistent and for the 
rest, it is not. Detailed analysis will show that the former is completely consistent. In particular these calculations 
correspond to the true minima of the free energy function. By contrast, the latter (non-self-consistent) portion 
is a consequence of a discontinuous transition which has "already" occurred at prior values of the parameters. In 
particular, the perturbative analysis is highlighting a local extremum and not the true minimum. 

We conclude this subsection with the derivation of the LSP-curve - as well as the introduction of notation 



that will be used throughout the reminder of the analysis. Assuming n = 0, Eq.(ll) is trivially satisfied and 



Eq.(lO) defines the "ambient" value of h which we denote by 6^: 



bB. := 



2e" 
1 + 2e" ■ 



Note that (6 = 6^, n = G = 0) is always a solution to the mean-field system. For simplicity we consider 
and /It as the relevant parameters for our system for the remainder of this paper. Let us now consider slight 
perturbations of h about and of n about zero. We thus write h = + A) with A ^ b^ and \n\ > with 



n <C 1 and obtain the following approximations by expanding Eq.(lO) to lowest order 



2A- ''^ 



1 



while Eq.(ll), written to a higher approximation than will be immediately necessary, gives us 



n nA 
bR bn 



Shu' 



(12) 



(13) 



We pause to observe that Eq.(13) and, in general, Eq.(ll), have the symmetry property that with all other 



quantities fixed, if n is a solution then so is —n. Thus, we might as well assume that n > 0. Indeed, we shall 
adhere to this convention throughout. Assuming now that /i is variable while ba is fixed, the n — ^ limit of 



Eq.(13l and Eq.(12| yields the tentative phase boundary 



1 

bR' 



(14) 



This defines the LSP-curve; the correspondingly tentative conclusion is that n > and b> bji occurs for n > fj,s 
while for fj. < i^s, n = b — bji = 0. However, the viability of these tentative conclusions depends, in a definitive 



fashion, on the value of fo/j. In particular an analysis of the higher order terms in Eq.(13) testifies that this 



picture cannot possibly be correct for 6ij < |; this is the subject of our next subsection. However, a more 



difficult analysis shows that this picture is indeed correct for 6_r > | 
subsection. First, we must attend to some necessary details. 



which is the subject matter of the final 



4-2.1. Preliminary analysis 

In this subsection we will establish some basic properties of the model such as the existence of high- and 
low-temperaturc phases along with various monotonicity properties. In particular we show that at fixed bn, the 
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quantity n, assumed to be non-negative, is non-decreasing with fi and strictly increasing whenever it is non- 
zero. For the benefit of our physics readership, such contentions might typically be assumed and consequently, 
the entire subsection could be skipped on a preliminary reading. However, it is remarked that in the normal 
(physics) course of events, such questions are most often settled by direct perturbative calculation. Even for 
continuous transitions, on some occasions, additional justification is actually required. Sometimes, as in the 
present work, when the transition is discontinuous, the relevant calculations simply cannot be done analytically 
and then, indeed, one must rely more heavily on abstract methods. 

In what follows, we shall work with the free energy function given by Eq.([6]) with G eliminated in favor of n 
according to Eq.([7]) and working with the parameters fi and For simplicity, this will be denoted by ^bit,iiib, n) 
but with subscripts omitted unless absolutely necessary. Thus $(6, n) is now notation for the function 

^b,i,t,{b,n) -.^ ~^nn^ - a{bR)b 

It is clear that the minimum of $(6, n) corresponds to the minimum of the original three variable free energy 
function $(n, 6, G). In the following, will use the notation n(/z) (with n(/i) > 0) as though this defines an 
unambiguous function. Of course in the case of phase coexistence, this will not be true. In general, then, n(/Lt) 
will stand for a representative from the set of minimizcrs at parameter value /i and all of the results in this 
subsection hold. We start with some elementary properties of the phase diagram generated by the corresponding 
minimization problem. 

Proposition 4.1 Consider $(6, n) with bji fixed and fi ranging in [0,oo). Then for all fi sufficiently large, 
$(6,n) is minimized by a non-zero n and for all /i sufficiently small, $ is minimized by {bj^,0). 

Proof. We begin with the assertion, gleaned from Eq.Q, that along the curves n = 0, $ is minimized by 
b = bjf. Thus we may pick any fixed, nontrivial no, with < no < b^, and it is sufficient to establish that 
^{fT-OTbif) < ^{0,bii) once fi is sufficiently large. However, the desired inequality is manifest for large fi since 
the only /x dependence in $ is in the term —^iiuq which is, eventually, in excess of the differences between the 
H independent term and $(0, 6/j). The second statement is proved as follows: since 6 = 0- which necessarily 



implies n = - does not minimize the free energy function, we may use the variable 9 :~ n/b so that Eq.(10l 
now reads 

&.^^^log(i± 

As is well known from the analysis of mean-field Ising systems (and can be established, e.g., by further differen- 
tiation) the above equation has only the trivial solution if 6/x < 1. Since b cannot be greater than one, the second 
statement has been proved - in fact whenever /i < 1. □ 

The above result establishes, in a limited sense, the existence of a phase transition. Here we will sharpen this 
result by proving that along the lines of fixed 6^, there is a single transition from n = to n > 0. This is an 
immediate corollary to the following lemma which we state separately for future purposes. 

Lemma 4.2 Let $^(6, n) denote the free energy function with bj^ fixed and fi (displayed) in [0,oo). Then the 
minimizing n{fi) , if unique, is a non-decreasing function of fi. More generally, if at various values of ^, has 
a minimizing set of n 's then, if fj,' > fi, the minimum of the minimizers at fi' is greater than or equal to the 
maximum of the minimizers at fi. Thus, in general any possible "choice" of n(^) is non- decreasing. 
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Proof. Let /i, /i' S [0,oo) with fi' > fi and let us denote by (b',n') a minimizing pair for and similarly 
for {b,n) at /i. The key observation is the meager /i-dependence of the fmrction Indeed, $^'(2:, y) — 

^ij-i^^y) — ~ ^J■)y^■ We do this twice: 

V(fe',n') = <i>^(6',n')-^(M'-M)[^f 
>$^(6,n)-i(M'-M)K]' 

= n) - - A^')"' - ^(a*' - (16) 

leading to ^^>{b', n') > ^^'{b, n) + — ^)(n^ — [n']^). This necessarily implies that [n'Y > since otherwise, 
the previous inequality would be strict implying that (6, n) would have been a "better minimizer" for than 

Using this result we may now show the following 

Corollary 4.3 Consider the mean-field model defined by the free energy function given in Eq.^. Then for 
each fixed bn £ (0, 1], there is a transitional value of fi, denoted by fixibR.), such that n = for fi < fix o,nd n > 
for > ^T- 



Proof. This follows immediately from Proposition |4.1| and Lemma |4.2| above. □ 
Also of interest is the following: 

Corollary 4.4 Consider the mean-field model defined by the free energy function given in Eq. Let 
denote any non-negative function corresponding to a minimizing n at parameter value fj, ( usually uniquely deter- 
mined). Then for ij,> ht, the function n{ij) is strictly increasing. 

Proof. It is seen that if u^jLt) = then the statement of this corollary is self-evident at ^ — jir- For the rest of 
this proof, we may simply assume that /i is such that n{fi) > 0. Suppose then that fi' > 11 and that n = n(fj,) is 
part of the minimizing pair (b{fi), n{ii)) at parameter value ^. Suppose further that at /z' the same n is also part 
of a minimizing pair. Then we claim that the b{p) is not the partner at /i' since given n' - purportedly equal to 
n - then b' is uniquely determined by Eq.(ll). Upon performing some algebraic manipulations the latter reads 
b' = n'/ tanh /x'n'. Thus, the equality n = n' would lead to 

6'^ ^ " ^ " =b 

tanh/i'n' tanh/i'n tanh /^n 

so that explicitly (6, n) cannot be a minimizer at parameter value fi' . Using the appropriate b' ^ b, we would 
have 

in contradiction with the fact that {b,n) is not minimized for the parameter value /z'. □ 
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4-3. A discontinuous Transition for bj^ < | 

The dividing point of bp = | along the LSP-curve n = l/bn is apparent from the higher order terms in 



Eq.(13). Indeed, supposing fj, — l/bn + e we obtain, with the additional aid of Eq.(12), 



1 



For bn > ^, Eq.(17) is consistent (and, as it turns out correct) but in the case of 6^ < ^ this equation alone 



precludes the possibility of a continuous transition. Indeed since we cannot have < 0, the only logical 



consequence of Eq. (171 is n = for fi > A*5(^_r), i-e., the transition occurs later. But the lower order term 
insisted that fj, = ^s(^-r) was the only viable candidate for a continuous transition. Thus: the transition cannot 
be continuous and, at least for 6/j < |, the preliminary assumption that n goes to zero continuously can no longer 
be sustained. In particular, for < |, perturbative analysis will never be valid because the relevant quantities 
will never be small. 

This leaves open the possibility of a transition at some ^T(^fl) that is different than /j.g. We shall show that 
< /^s as a direct consequence of the following: 

Proposition 4.5 Consider the mean-field model defined via the free energy function given by Eq.^. Then, if 
bji < 7^ at jjL = HsibB.) the quantity n is strictly positive. 

Proof. We expand the free energy function <I>(n, 6) - with G eliminated via Eq.[7]- about b — bp;, n = and 
along the curve = /i5(&fl). The convenient variables are now chosen as 6 = 6_r,(1 + A) and n = bum. We first 
note that all odd terms in m must vanish. In addition, the term linear in A vanishes due to the stationarity of $ 
along the curve /i = ^J■sibR) and, as it turns out, so does the term which is quadratic in m. This leaves us with 



^b,n) = ^bR,0) + hR 



-to" H — A^ - m^A 

6 1 - Ob 



+ . 



Examining the quadratic form in the variables to^ and A, the condition for a local minimum is that 

1 1 1 

61 -bn ^ 4' 

i.e., bji > ^. We return to bu > ^ in the next subsection. Of current relevance is the fact that for bp < |, the 
curve fj, = b~j^ is of a saddle point nature. This implies that there is a direction of decrease which, as is easily 
seen, is optimized, in the physical direction, when = 3A. It is concluded that under the stated conditions, 
we can produce a pair (6, n) with > (and b > bp) such that the free energy for the non-trivial pair is lower; 
we just make the corresponding objects small enough to withstand the higher order corrections. Thus the actual 
minimum also must occur for non-trivial values of the n observable. □ 

We now have 

Theorem 4.1. Consider the mean-field system defined by the free energy function as given in Eq.(|6|. Then for 
bn. < gi there is a discontinuous transition at some positive nxibR) < Ais(^fl)- 



Proof. That a transition occurs at some fix > is the statement of Corollary |4.3[ Moreover, Lemma 4.2 and 
the above analysis implies < Ms- The discussion prior to Proposition |4.5| demonstrates that at jj, — fiT, the 
quantity n is already positive. It only remains to show that the inequality relating fis and /it is strict. To 
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this end, let us reimplement the heretofore unnecessary notation for the full dependence of the free energies on 
parameters. We have learned that for bp < |, there is an n^, > and a 6^, (with bi, > bji) such that 

FMF{bR,^J.s) = $b„,ps(6*,n*) < $b^,,,g(6_R,0). 

Invoking Lemma |4.2[ it is now sufficient to show that there is a > such that for some nonzero n, and some 
b, the inequality n) < <&6j^,^s-5^(fo_R, 0) can be shown to hold. Once again, the key is the simple 

dependence of the free energy functions on the parameter /i. Indeed, using rii, and 6* as trials, we obtain 

while 0) = $f,^_^g 0) < F^^^^^. Thus, the desired inequality will indeed hold for all 6fj, 

sufficiently small. □ 

4.4- A continuous transition for bji > ^ 

The starting point in our analysis is to show that at the purported critical curve, the quantity n actually 
vanishes. 

Proposition 4.6 For bn > ^ and /i — 6^^ =: /ig, the unique solution to the mean-field equations is n = with 
b = bR. In particular, $(,^^^3 (6_r, 0) < $6^,^3(6,71) for any {b,n) ^ (6^,0). 

Proof. Assuming n > the agent fraction b can be eliminated in favor of the ratio 

Note that while this is the same substitution as before, here it is b rather than n that is being eliminated. 
Notwithstanding, 9 still satisfies < 6* < 1. In these variables, the mean-field equations, Eq.(lO) and Eq.(ll| 
respectively become 

RO 

n= , 18 

n = 6_R Arctanh 6* (19) 
where in the above, R := 6_r/(1 — b]^). Let us now define £{9) as 

e(0) 1. = (1 + ^ + 

bB. R + B + R + Q ' 

where Q = Q{9) := \/l — 6'^. To prove the current proposition we need to show that for all 6 > 0, 

Arctanhe > £{9) 

demonstrating that there cannot be a non-trivial solution to the mean-field equations under the conditions 
stated. Note that for < ^ 1 the desired inequality can be explicitly demonstrated. In general, it is sufficient 
to show, for < 6* < 1, that £'(9) < 1/(1 - 9'^), i.e., in the Q variable that 

^ {1 + R){R + Q + {1-Q^)/Q) 
g2 ^ (R + Qy 
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Although both sides diverge as Q — ^ the divergence on the left hand side is clearly stronger so we actually 
only need consider Q > limiting us to Q € (0, 1). After some manipulation, the inequality we need to prove is 
equivalent to 

{R + Qf > (1 + i?)(i?Q2 + Q) = (1 + R){RQ-^ + g2) + {l + R){Q- q2). 
That is, we now wish to show 

R{R + 2Q + RQ){1 - Q) > (1 + R)Q{1 - Q). 
Since Q ^ 1, the above is equivalent to 

R^ + QR + R^Q > Q. 

Finally, since also Q < 1 it is enough to show that 2R^Q^ + RQ > Q i.e., that 2R^ + R > 1 which occurs for 
i? > i. This corresponds to > i. □ 

We can finally show 

Theorem 4.2. Consider the mean-field Gl-system defined by the free energy function given in Eq.(|6]). Then, 
for bji > as a function of with bu fixed, there is a continuous transition at ^ = /is = I.e., n(/Lt) = for 

fj, < and n(/i) > for the n-component of any minimizing pair (n(/x), &(//)) while, if /x 4- fJ-s, it is found that 
i 0. 

Proof. We will marshal the facts at our disposal and then proceed in a more abstract vein than has been the 
case in the more recent of our arguments. In what is to follow, and the corresponding 6(/i) is, once again, 
notation for a minimizing pair without any claims to uniqueness. By the preceding proposition, we know that at 



fi = fis, the quantity n{fi) is unambiguous and vanishes for fi < fis hy Lemma 4.2 Conversely, for fi > fis 'Vfe 



may write, adhering to the notation in the proof of Theorem 4.1 our usual expression: 



1. 



For oc b — bfj 1 from Proposition 4.5 we know that the quantity ^bR,fj.s (^i agrees with ^bR,tis i^R^ ^) 
quartic order in n. Thus allowing <C 1 with n'^{^ — fig) 3> n*, {b — fe/j)^ we find a non-zero n corresponding to 



a free energy lower than that of ^bR.p,s i^R^ 0)- Therefore, again by Lemma 4.2 we have n(/i) > for all ^ > fj,s- 
It remains to establish that n J, as J, /ig. Note, that along any decreasing sequence of /i's the corresponding 
possible n's must be monotone by Corollary |4.4| - or even Lemma [4. 2| - and hence n I as ^ I fis- Now let us 
suppose otherwise: that for some sequence of /i's decreasing to fis there is an associated sequence of minimizers, 
(6(/i), n(/i)) that has Ti{fi) ln^, > 0. Let b^, denote the associated limit for the b{fi) along a further subsequence 
if necessary. Since 

$f,„,^(6(/i),n(/i)) < $6^,;.(5fl,0) EE $6„,^s(5fl,0) 

we would have, by continuity, ^ba,ns — ^ba,p,s {bR, 0) indicating that at fi = fis, there is a minimizer with 

positive magnetization in contradiction with Proposition |4.6| above. It follows that, under the stated condition 
^, the limit of n{fi) is zero as fi fis while it vanishes below and is positive above. By this (and any other) 
criterion, the transition at fis is continuous. This completes the proof. □ 
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5. Discussion 

In this work, we have formulated a lattice model for gang territoriality where red and blue gang agents interact 
solely through graffiti markings. Using a contour argument, we showed that a phase transition occurs between 
a well mixed, "high-temperature" phase and an ordered, "low-temperature" one as the coupling parameter J 
between gang members and graffiti becomes stronger while the graffiti evaporation parameter A decreases. In the 
mean field limit of all-to-all lattice site couplings, we can also identify the tricritical point in phase space that 
distinguishes the occurrence of a continuous phase transition from a first order one. We find this point to be 
located at b^. = 1/3 which corresponds, in terms of the original variables of the problem, to the gang proclivity 
term a = — 21og2. In particular, for > 1/3 the phase transition is continuous and occurs at ii — l/bn.. Thus, 
in the mean-field limit, for fixed a > —2 log 2 the ordered "low temperature" phase arises for > A/(e~" + 2), 
and the "high temperature" one is attained on the other side of this inequality. The transition between the 
two occurs in a continuous manner across the = A/(e~" -I- 2) locus. In the opposite case of bji < 1/3 (or 
a < —2 log 2) the phase transition is discontinuous. Here, we also are able to prove that the transition between 
high and low temperature phases occurs not at /i = l/bji, but rather along the fj, = fj,T < ^/bR curve, so that the 
phase change occurs earlier in J and along a separatrix = < A/(e~" -|- 2). 

In the context of gang-graffiti interactions, we may identify the low temperature, clustered phase as pertaining 
to a high level of antagonism between between rival gangs, where segregation leads to conflict along boundaries. 
Vice versa, the high temperature, well mixed configuration can be interpreted as a peaceful state, where despite 
different affiliations, gang members share the same turf. Our mean field results indicate that the confrontational 
state is surely attained, whether in a continuous or first order manner, for > A/(e^"-f 2), which represents high 
gang-graffiti territoriality J, low external intervention in graffiti removal A and high proclivity a for individuals 
to become gang members. Gang clustering can be avoided by intervening in all three directions: by externally 
eliminating graffiti (A), but also, from a deeper sociological point of view, by decreasing the lure of graffiti tags 
or of joining gangs in the first place (J, a). The emergence of a (continuous or discontinuous) phase transition 
shows that it is possible to obtain segregation in a lattice model without invoking direct agent-to-agent coupling; 
it is certain that adding such coupling terms to the Hamiltonian would allow for even more favorable segregation 
conditions. 

Although our work was conceived within the context of gang interactions, the proposed model Hamiltonian 
and the tools used are general enough that our fundamental results may be applicable to several other contexts 
where territoriality is played out through markings and not through direct contact between players. Many 
animals, among which wolves, foxes and coyotes, are known to scent -mark their territories as a way of warning 
intruders of their presence and to exchange internal communication [41 1. At times, buffer zones can originate 
between distinct animal clusters where prey species, such as deer or moose, may thrive |63| . Insects, such as 
beetles and bees, are also known to avoid previously marked locations as a way to optimize foraging patterns. 
Similarly to the role of gang graffiti markings, foreign scents lead "others" to retreat from already occupied turf 
or visited patches. Our work also applies to these contexts. Although some stochastic treatments have been 
recently presented [33] , classical ecological studies of territoriality are usually carried out via reaction-diffusion 
equations where focal points such as dens, burrows or nests are often included [321 ISO], leading to segregation. 
Within this work on the other hand - whether first order or continuous - agent clustering is a natural consequence 
of a probabilistic treatment without the need to include any anchoring sites. Finally, we are able to connect local 
microscopic parameters - J, K, X,a - to the emergence of large scale territorial patterns, be they gang clusters 
or animal groupings. 
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Appendix A. 

Here we present a brief description of the informal simulations mentioned in Section|4] We consider a 100 x 100 
square lattice with periodic boundary conditions initialized at i = so that each site is populated with either 
red or blue agents, or a mixture of both. We assume a random distribution of 10^ blue and 10^ red agents and 
do not impose any restriction on the number of individuals on each site, so that multiple agents can occupy the 
same location at any given time. Initial conditions are completed by assuming that at i = there is no graffiti 
present. At each time step of the simulation agents leave their graffiti on-site with a probability pm that depends 
on the current graffiti level. In particular, = 0.1 if the site is not marked by the opposite gang, and Pm = 1 
otherwise. The agent then moves to any of its nearest neighbor sites j with probability e"^^ /^^ e^^ , where gj 
is the amount of the opposite gang's graffiti at location j. Agents will thus preferentially relocate to nearest 
neighbor sites tagged by the least amount of the opposite gang's graffiti. Finally, at each site, graffiti is removed 
according to a probability Pg. Similarly to the number of agents, we impose no restriction on the amount of 
graffiti at each site. 

While the simulation rules described here are similar in spirit to the model we analyze in this work, they do 
not directly lead to the Hamiltonian in Eq.[T] These informal simulations however provided us with a playing 
ground to investigate any phase transitions that may take place upon varying relevant parameters, such as the 
graffiti removal probability pg. For example, in Table A.l[ we track the dynamic progression of two sets of 



parameters. In the top row we set pg — 0.25 so that 75% of the graffiti is retained at each iteration, while in the 



lower one we set = 0.75 so that only 25% of the graffiti is kept. In Table A.l the left and right hand side plots 
show agent and graffiti distributions, respectively. Red and blue pixels indicate site occupied by respective gang 
agents, black pixels represents no agents, and magenta shades indicate coexistence of both red and blue agents. 
Just as in the main body of this paper, we do not include any direct coupling between red and blue agents who 
interact only via the graffiti field. Similarly to what wc later found in the main analysis, the degree of persistence 
of the graffiti field - which can be related to A in the Hamiltonian in Eq.[l]- yields different qualitative behaviors 
and, if sufficiently large, may lead to aggregation patterns with distinct red and blue phases. The emergence of 
separate clusters from these simulations motivated the more extensive study presented in this work. 
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Table A.l; Snapshots of a Monte Carlo simulation of gang dynamics on a 100 X 100 square lattice with periodic boundary conditions. 
For each image in the sequence, the upper left panel represents gang agent populations while the upper right panel is the corresponding 
graffiti distribution. Iteration time is measured in arbitrary units. At t = 0, 10^ red and 10^ blue gang members are placed at random 
on the lattice with possible overlaps. Magenta indicates a mixture and black indicates a void in gang agents or graffiti. Agents tag 
their sites with probability pm =0.1 if the site is not marked by the opposite gang's graffiti and with probability pm = 1 otherwise. 
In the upper panels of the table entries, pg = 0.25 so that at each time step, graffiti will persist with a 75% possibility. The lower 
panels, where pg = 0.75, mirror the upper ones but with a much lower graffiti persistence, of 25%. Note the different outcomes of the 
simulations at long times; when graffiti is allowed to persist longer, segregation occurs with the formation of islands of red and blue 
gangs. In this work, just as in our current model, there is no direct interaction between gang members, underlying the importance 
of the graffiti field as an indirect coupling between agents. 
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